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ABSTRACT 

We have performed multivariate statistical analyses of photometric and chemical abun¬ 
dance parameters of three large samples of stars in the globular cluster ui Centauri. 
The statistical analysis of a sample of 735 stars based on seven chemical abundances 
with the method of Maximum Parsimony (cladistics) yields the most promising results: 
seven groups are found, distributed along three branches with distinct chemical, spatial 
and kinematical properties. A progressive chemical evolution can be traced from one 
group to the next, but also within groups, suggestive of an inhomogeneous chemical 
enrichment of the initial interstellar matter. The adjustment of stellar evolution mod¬ 
els shows that the groups with metallicities [Fe/H]> -1.5 are Helium-enriched, thus 
presumably of second generation. The spatial concentration of the groups increases 
with chemical evolution, except for two groups, which stand out in their other prop¬ 
erties as well. The amplitude of rotation decreases with chemical evolution, except 
for two of the three metal-rich groups, which rotate fastest, as predicted by recent 
hydrodynamical simulations. The properties of the groups are interpreted in terms of 
star formation in gas clouds of different origins. In conclusion, our multivariate anal¬ 
ysis has shown that metallicity alone cannot segregate the different populations of ui 
Centauri. 

Key words: Globular clusters: individual: ui Centauri (NGC 5139) - stars: evolution 
- stars: AGB and post-AGB - methods: statistical 


1 INTRODUCTION 


It is now well established that all globular clusters (or at 
least all those surveyed in depth) harbour multiple popu¬ 
lations of stars. Evidence for this has been found in their 
spectra, which reveal different abundances of elements, and 
in their color-magnitude diagrams (hereafter CMD), where 
distinct main sequences or subgiant branches are visible. 
While the exact mechanism is still debated, the generally 
accepted explanation is that these distinct populations are 
successive generations of stars, the later ones having formed 
from the ejecta of the earlier generations. Other explanations 
for these multiple populations, such as mergers of globular 
clusters, original chemical inho mog eneities or p ollution by 
interacting massive binaries (Bastian_etaL| 2013 1. are how¬ 
ever not excluded (see e.g. iGratton et all 20121 . for a recent 
review). 

ui Centauri, the most massive globular cluster of our 
Galaxy, stands out in this respect. While a modest spread 
in [Fe/H] (0.2 dex or less) has been found in some of the most 
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massive Galactic globular clusters, t h e spread in uj Cent auri 
is more than 1.5 dex (see [Johnson fe Pilachowski]|2010l . for 
a recent update), leading to speculations that it might be 
the remains of a tidally stripped dwarf galaxy. Being only 
5.2 kpc from the Sun, w Centauri has been extensively in¬ 
vestigated and the properties of its stellar populations well 
documented. 

Most studies use the CMD to distinguish the dif¬ 
ferent stellar populations of globular clusters, and gener¬ 
ally identify stellar sequences with slightly different ages, 
met alli cities and more rarely different helium abundances. 
iBedin et al.l (120041 ) identify at least two main sequences in 
the CMD of u> Centauri and offer several possible expla¬ 
nations in terms oj) c hemic al abun dances or pollution by 
background stars. ijoo &: Leel (120131 ) identify five different 
populations of increasing metallicity and helium abundance 
in the CMD of u> Centauri: a metal-poor one (13.1 Gyr), 
three of intermediate metallicity (13.0, 12.0 and 11.4 Gyr) 
and a metal-rich one (11. 4 Gyr). Focu sing on the subgiant 
branch of uj Centauri, IVillanova et alj (120141 ) find 6 popu¬ 
lations with internal age spreads of over 2 Gyr, and some 
stars that are merely 4 or 5 Gyr old. 

Another way of separating different stellar populations 
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is by their chemical abundances, using spectroscopy. Since 
this is expensive in terms of teles cope ti m e, few stu dies 
are based on large samples of stars. 1 Johnson fc Pilachowskil 
(|2010f) , whose results are used in the present paper, de¬ 
termined abundances of various elements of 855 stars on 
the red-giant branch of uj Centauri from high-resolution 
(R=20000) spectroscopy. They fit gaussians to the metallic- 
ity distribution of these stars, and distinguished four popula¬ 
tions in this way: RGB-MP, RGB-Intl, RGB-Int2+3, RGB-a 
(metal-rich). 

In standard mono-metallic clusters, different genera¬ 
tions of stars can also be segregated via their location in 
the Na- O plane. Using chem ical abundances in 19 globular 
clusters ICarretta et al.l (|2009l J derived an empirical criterion 
for separating the first (or primordial) and second genera¬ 
tion of stars. The former, which contains about one third of 
the stellar populations, is limited to [Na/Fe]< Na m i„ + 0.3. 
In the latter they distinguished an intermediate component, 
which contains about 60 to 70% of the stellar populations 
and an infrequent extreme one, which is mainly present in 
more massive clusters. As will be discussed below, the mix¬ 
ture of stellar populations in uj Centauri is too complex for 
their criterion to be applicable. 

There is thus ample evidence for multiple populations 
in uj Centauri, as well as in other globular clusters, but the 
methods for distinguishing them are rather empirical, us¬ 
ing generally o ne parameter ( the metallicity) or two (those 
of the CMD). IGratton et all (1201 lh were the first to ap¬ 
ply a method of multivariate analysis to separate the stel¬ 
lar^ popula tions o f globula r clu sters. They used the data of 
I Johnson fc Pilachowskil ( 201f)f b and applied a k-means al¬ 
gorithm to four parameters of 797 stars: [Fe/H], [Na/O], 
the average of the abundances of Si, Ca, and Ti, and the 
abundance of La. They opted for 6 groups (in k-means the 
number of groups is an input parameter), which they further 
rearranged into three main groups according to metallicity. 
They did not consider these three main groups as successive 
generations, because the Na-O anticorrelation is present in 
the first two, suggesting the presence of first- and second- 
generation stars in both. Na and O are correlated in the 
third main group, unlike in any other globular cluster, but 
the authors found no evidence for an external origin of these 
stars._ _ 

ISimpson et al.l (|2012| J used low-resolution (R=1600) 
spectroscopic data of 221 red giants in uj Centauri 
to derive their effective temperature, gravity and four 
ab undances, to w hich they added two abundances from 
[Johnson fc Pilachowskil (l2010ll . They performed a k-means 
analysis on six parameters ([Fe/H] and the abundances of 
C, N, O, Na and Ba) under the assumption of four groups. 
The analysis revealed a low-metallicity group (which they 
identify as the primordial population), two intermediate- 
metallicity ones differing in their O and Na abundances 
(one of them presumably helium-enriched) and a high- 
metallicity group (in which Na and O are correlated). 
The two intermediate-metallicity groups are blended in the 
CMD; when combined, they show the well-known Na-O anti¬ 
correlation. The authors identified their low- and high metal¬ 
licity groups with the corresponding ones of IGratton et alj 
(1201 ill . In one of their proposed scenarios the intermediate- 
metallicity group of high O abundance would be the last 
one to form, from the gas of the primordial generation, af¬ 


ter a gas-sweeping passa ge through the Gala ctic plane. In a 
subsequent study dSimpson fc Cottrelll l2013| j they classified 
848 stars, 557 of which are new ones, into these four groups 
through a kind of decision tree. They concluded that forma¬ 
tion models of globular clusters and their stellar populations 
do not provide a consistent picture for the four groups. 

In the p re sent pa per we repeat the analy ses of 
Sim| 3 Son_etjiL| d2012h . ISimpson fc Cottrell (l2013li and 
Gratton et al.l (1201 il l using a different method of multi- 
variate analysis, that of Maximum Parsimony (cladistics, 
iFraix-Burnet et al.ll2006bl la). which requires no a priori as¬ 
sumption about the number of groups, and which does pre¬ 
cisely what the other statistical methods are unable to do: 
order the groups chronologically. In all this paper, the cladis- 
t ic a nalyses were performed with PAUP*4.0bl0 (jSwoffordl 
l2003ll with the parameters discretized into 30 equal-width 
bins, which play the role of discrete evolutionary states. 
More details are given in Appendix [A] We also apply the 
k-means algorithm to the data, in order to test the repro¬ 
ducibility of statistical analyses with k-means, and compare 
the results give n by the two methods (se e more details on 
the methods in IFraix-Burnet et al .1120121) . 


2 STATISTICAL ANALYSES 


2.1 The Simpson et al (2012, 2013) samples 

We first analyse the sample of 221 stars of ISimpson et all 
(|2012|) , using the same six abundance parameters as them. 
In partitioning analyses like k-means, the number of groups 
must be given as an input. To estimate the optimum num¬ 
ber of groups, we have used several statistical criteria: the 
Calinski-Harabasz criterion and the Simple Structure Index 
through the function cascadeKM of the R package vegan 
(|R Core Teaml 1 2 01-ll: Oksanen et al. 2013l b and the jump 
technique from ISuear fe Jamesl (120031 ). Even though these 
tests are mainly indicative, we find the optimal number to 
be two, or nine for one of the tests, but n ever fo ur. As a 
rationale behind the choice of four groups, ISimpson et ali 
d2012li state that their data are not good enough for divid¬ 
ing their sample into more groups. If we impose four groups, 
the results of our k-means analysis disagrees with theirs in 
the number of stars in each group: 9, 52, 60 an d 100 in our 
case, and 24, 49, 62 and 86 in theirs. ISimpson et al.l d2012h 
do not provide much details on how they run the k-means 
algorithm, and it is possible that their study suffers from 
one well known caveat of the k-means: the result depends a 
lot on the initial seed. To get rid of this depen dence, we used 
the function kmeans of the package stats in R dR Core Teaml 
l2014lh which repeats the analysis many times (we chose 1000 
repeats) with a different seed and selects the best cluster¬ 
ing. No further comparison is possible because the details of 
their groupings are not provided. 

Our cladistic analysis on the sample of 221 stars and 
six abundances gives six main groups which are not in good 
agreement with the four groups of our k-means clustering 
described above, but in much better agreement with a k- 
means analysis with six groups. 

We again do not find any statistical evidence in favour 
of an optimal number o f groups larger than two w ith the 
sample of 848 stars from ISimpson fc Cottrelj d2013l b In ad¬ 
dition, they do not perform any clustering analysis. They 
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Table 1. The 24 parameters and the number of undocumented values from the compilation by iGratton et al.l 12011 *). 


Parameter 

Unknowns 

Parameter 

Unknowns 

V magnitude 

0 

[Si/Fe] 

34 

B - V 

0 

[Ca/Fe] 

1 

J magnitude 

0 

[Ti/Fe] 

32 

H magnitude 

0 

[a/Fe]= mean ([Si/Fe], [Ca/Fe], [Ti/Fe]) 

0 

K magnitude 

0 

[Na/O] 

14 

T e fr 

0 

[Scl/Fe] 

441 

l°g(g) 

0 

[ScII/Fe] 

100 

[Fel/H] 

0 

[Sc/Fe]= mean([ScI/Fe], [ScII/Fe]) 

37 

[FeII/H] 

164 

LaII/Fe 

48 

[Fe/H]= mean([FeI/H], [FeII/H]) 

0 

[Ni/Fe] 

52 

[O/Fe] 

7 

[EuII/Fe] 

721 

[Na/Fe] 

10 

[Al/Fe] 

723 


classify by hand 557 stars not present in the first sample 
by comparing the four parameters they found to character¬ 
ize best the four groups from their first study. This method 
has serious limitations, because the larger sample covers a 
broader domain in the CMD. 

2.2 The Johnson & Pilachowski sample: four 
parameters, 797 stars 

1 Johnson fe Pilachowski! (feoioh provide eleven abundance 
parameters for their sample of 855 stars, to which 
IGratton et al.l (|201ll ~) added photometry, T e fr and log(g) 
leading a total of 24 parameters listed in Table |T] If we con¬ 
sider that some of them ([Fe/H], [Sc/Fe], [a/Fe]) should be 
replaced by mean values to compensate for undocumented 
values and two ([Al/Fe] and [EuII/Fe] ) should be discarded 
for lack of data, we end up with 15 parameters usable for a 
clustering analysis. 

To complete the table, we computed the global metal- 
licity Z using log Z = [Fe/H]- 1.72125 from ICarraro et al.1 
d 19991'). The range of Z is 0.0004 to 0.003. 

iGratton et al.l (|201lh selected four parameters based on 
a priori physical arguments for detecting the different gen¬ 
erations: 

• [Fe/H] is assumed to be representative of the overall 
metallicity, 

• [Na/O], which is representative of the location of stars 
along the Na-0 anticorrelation, and is likely to be correlated 
with the He abundance, although this relation is presumably 
not linear, 

• [a/Fe] is the average of Si, Ca, and Ti, all mainly pro¬ 
duced by a-capture reactions in massive stars, later explod¬ 
ing as core-collapse SN, 

• The abundance of La ([LaII/Fe]), which is an n-capture 
element that in the Sun is mainly produced by the s-process. 

After removal of undocumented stars for these four 
parameters, t he 797-star dataset includes the 221 of 
ISimpson et al.l (120121 ') but [Fe/H] is the only parameter in 
common between the two k-means studies. 

We perf o rmed the same k-means analysis as 
IGratton et al.l d201ll ~) and found exactly the same dis¬ 
tribution in size of the six groups. Unfortunately, the 
statistical tests prefer two groups, maybe four, but not six. 
There is thus no statistical justification for the six groups, 


which means that the algorithm is f orced to art ificially 
divide the sample into six groups. IGratton et alj d201lh 
found that if they choose three groups for the k-means 
analysis, the differences between the groups are largely 
dominated by [Na/O]. Since they do not find this result so 
informative, they thus increased the number of groups to 
six. 

The cladistic analysis of the same data set with the same 
four parameters produces three main groups that include 
many small ones. The three groups are easily distinguishable 
in [Fe/H], [LaII/Fe] and [a/Fe], and one has a higher aver- 
age [Na/Q[. Note th at, in contrast to the three groups that 
IGratton et al.l d201 if) found by k-means clustering, the three 
groups resulting from our cladistic analysis are not discrim¬ 
inated by [Na/O]. Nevertheless, we do not find any really 
new in formation with respect to the study by IGratton et all 
(l201lh . We believe that this is due to the choice of parame¬ 
ters and of number of groups, a choice which is biased by a 
prioris coming from the physics of stellar evolution, so that 
the settings of the analysis are not optimal from a statisti¬ 
cal point of view. More precisely, [Na/O] appears in the first 
component of the Principal Component Anal ysis, it exp l ains 
the axis of greatest variance in this sample (IGratton et al.l 
l201ll) . This does not mean that it is the most discriminant 
one. 

To test the influence of [Na/O], we considered Na and 
O separately, so we use the five parameters [Fe/H], [Na/Fe], 
[O/Fe], [LaII/Fe] and [a/Fe] to find the optimum number 
of groups for a k-means analysis. The result is more satis¬ 
factory since the optimal number is clearly three or four, 
as compared to only two when one looks for the optimum 
grouping using the initial 4 parameters. This trend is even 
stronger if we add log(g). I n othe r word s, the choice of pa¬ 
rameters made by IGratton et al.l (1201 ll*) was probably not 
the best one for the task. We thus decided to perform a 
fully objective analysis by selecting the parameters from an 
entirely statistical point of view. 

2.3 The Johnson & Pilachowski sample: seven 
parameters, 735 stars 

We now repeat th e clu stering analyses on the 
Ijohnson fc Pilachowskil (|2010| J sample with an optimised 
choice among the 15 usable parameters. In a preliminary 
step, we checked all the correlations, using all the scatter 
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plots together with a detailed examination of the results 
from a Principal Component Analysis, looking carefully for 
the influential parameters within each component. 

After this first step, we kept the V magnitude and dis¬ 
carded H, J and K because the correlations among them are 
strong and redundancy reasonably present. Also B — V is 
strongly correlated with T e g and very probably redundant, 
so that we discard it as well. Following the result of Sect. 12.21 
we did not include the ratio [Na/O]. 

We are left with ten parameters: V, log(g), T e ff, [Fe/H], 
[O/Fe], [Na/Fe], [Sc/Fe], [Ni/Fe], [Lall/Fe], [a/Fe], How¬ 
ever, V, log(g) and T e fi are very much correlated, especially 
the two first ones, and largely dominate the first Principal 
Component. In addition, a cladistic analysis using the ten 
parameters results in a very linear arrangement of the tree 
according to these three parameters. This indicates a clear 
redundancy. We thus must discard at least one of them. The 
two parameters Vmag and T e g have a complicated evolu¬ 
tionary behaviour in the considered region of the HR di¬ 
agram where the evoluti onary tracks (and t he iso chrones 
as well) show loops (e.g. lOirardi et al.ll200fl 120041 1 . These 
back and forths preclude these parameters from being reli¬ 
able indicators of a population or a generation. Note that 
this argument is especially valid for cladistics, but also for 
partitioning techniques if the goal is to gather objects with 
the same history. 

We finally retain the following seven abundance param¬ 
eters: [Fe/H], [O/Fe], [Na/Fe], [Sc/Fe], [Ni/Fe], [Lall/Fe], 
[a/Fe]. Some of these parameters have up to 52 undocu¬ 
mented values (Table 0, that is less than 7%. This is not a 
problem for cladistics, but k-means, like all distance-based 
approaches, cannot be used in such a situation. It is possi¬ 
ble to replace these missing data for instance by the mean 
of the each parameter, but we decided to consider a fully 
documented subset. The fina l sample has 7 3 5 sta rs, which 
is slightly fewer than used bv lOratton et al~. '. (12011 /1. 

The optimal number of groups for the k-means analysis 
is found to be either 2, 3 or large like 35 depending on the al¬ 
gorithm. For comparison with the cladistic result below, we 
set the number of groups at three or seven. The distribution 
of stars among the groups is given in Table [2] 

The cladistic analysis of the 735 stars with seven pa¬ 
rameters produces a tree that is quite robust (each node oc¬ 
curs in more than 88% of 9517 equally parsimonious trees). 
Three main branches, defining three main groups, are iden¬ 
tified on Fig. 0 as OC1, OC2 and OC3. The latter two can 
be divided into smaller groups: OC2a, OC2b, OC2c, OC3a, 
OC3b and OC3c, making seven groups altogether. The total 
number of objects in each group is given in Table [2] Many 
other sub-branches could define groups but their number of 
objects would be rather small, and the subsequent interpre¬ 
tation of the tree does not seem to justify this. Note that 
from a strictly cladistic point of view, the OC2a and OC3a 
groups may not represent a fully homologous ensemble of 
objects since they do not include all descendants from the 
uppermost node. They might be gathered with one of the 
two other groups (for instance OC2a and OC2b could make 
a single group from which OC2c diverges), but this cannot 
be decided without the analyses of the parameters done in 
the remainder of this paper. 

The contingency table (Table 0 provides a detailed 
comparison between the results of the k-means and cladistic 


OC2c 


t 



Figure 1. Unrooted tree obtained with cladistics on 735 galax¬ 
ies described by seven parameters (see Sect l2.3l l. There are three 
main groups (branches OC1, OC2 and OC3). Groups OC2 and 
OC3 are each divided into three groups labelled a, b and c (see 
text). The color codes of the groups will be used throughout the 
paper. 


results. For three groups, the k-means result is in fairly good 
agreement with the cladistic analysis since between 70 and 
95 % of the objects of a group are found in a same group 
of the other technique. For seven groups, the agreement is 
poorer. This is not surprising since seven is not found to 
be an optimal number of groups for the k-means analysis. 
Strictly speaking, cladistics may identify five groups (see 
above), but a k-means analysis with this number of groups 
does not improve the agreement. Maybe a higher number of 
groups could be pertinent, as suggested by the k-means tests 
of the optimal number of groups and in agreement with the 
numerous small subbranches that can be identified on the 
tree in Fig. 0 but the number of stars within each group 
would be too small for an interesting astrophysical interpre¬ 
tation. In addition, the agreement would certainly not be im¬ 
proved by merely subdividing already discrepant groups. We 
conclude that the agreement for three main groups is very 
much satisfactory, and that the k-means approach has some 
difficulties to find more subtle clustering in this dataset. 

The cladistics result presented in this Section (Fig. 0, 
using 735 stars and seven parameters, is analysed in the 
remainder of this paper. 

3 PROPERTIES OF THE SEVEN GROUPS 
3.1 Chemical evolution 

Since the cladistic analysis is based on abundance parame¬ 
ters only, one expects the groups to trace the chemical evo¬ 
lution of the stars in u Centauri. This is indeed what Fig. lBlI 
shows, a progressive chemical evolution, as [Fe/H], [Na/O], 
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Table 2. Contingency table comparing the k-means and the cladistic clusterings of the sample with 735 stars described by seven 
parameters (Sect. Ol i. The left part of the table considers the three-group hypothesis, and the right part considers the seven-group 
hypothesis. 



three groups 



3 


seven 

groups 







Cladistics 

group 

Number 
of objects 

k-means 

1 2 

Cladistics 

group 

Number 
of objects 

i 

2 

k- 

3 

means 

4 

5 

6 

7 

OCl 

230 

195 

33 

2 

OCl 

230 

135 

45 

3 

46 

0 

0 

1 

OC2 

265 

66 

190 

9 

OC2a 

119 

9 

53 

4 

38 

12 

0 

3 






OC2c 

99 

1 

31 

54 

13 

0 

0 

0 






OC2b 

47 

0 

1 

1 

45 

0 

0 

0 

OC3 

240 

18 

30 

192 

OC3a 

47 

8 

0 

9 

14 

10 

0 

6 






OC3c 

143 

1 

1 

0 

0 

61 

41 

39 






OC3b 

50 

0 

1 

6 

0 

3 

0 

40 


0C1 



-1.0 -0.5 0.0 0.5 

[O/Fe] 








Figure 2. [Na/Fe] vs [O/Fe] anticorrelation. It is present globally, 
and to some extent in OC2a and OC1. The correlation is in OC3c. 

[Al/Fe], [a/Fe] and [Lall/Fe] increase broadly from the OC1 
to the OC2 and OC3 main groups. 

There are distinctive features that explains the divi¬ 
sion of the three branches into seven groups (Figs. IB 1 1 and 
IB2I) . Within OC2, the groups OC2a and OC2b differ in 
[Ni/Fe], La/Eu and marginally in [a/Fe] and [Al/Fe]. In 
OC2b [Ni/Fe] increases strongly with rank. The group OC2c 
has lower T e fj and log(g), and higher [Fe/H] and [a/Fe] than 
the two others. It also has brighter K magnitudes. The group 
OC2c looks similar to OC3 in that [Fe/H] and [a/Fe] in¬ 
creases with rank, but it has a high [O/Fe]. In other words, 
the mechanism for producing the Na-0 anticorrelation is not 
at work in OC2c. 

The groups within OC3 have the same high [Na/Fe] and 
[a/Fe], but differ in some of their abundances, OC3a having 
lower [Fe/H], higher [O/Fe] and [Ni/Fe] than the others, and 
OC3c having the highest [Lall/Fe] of all. 

When examining Fig. Al, one should keep in mind that, 
cladistics being a most parsimonious approach, the ordering 


Figure 3. [Ni/Fe] vs [Fe/H] (weak) anticorrelation. It is present 
in OC3b and OC3c, and marginally in OC1. The groups OC2a 
and OC2b are clearly separated in this diagram. 

of the stars within branches favors a regular increase of the 
parameters used for the computation, if they are compatible. 
This explains why the dispersion is high for the parameters 
not used in the cladistic analysis (photometry, log(g), T e fi). 
Incidentally, in the main group OC1, we note a correlation 
between V, K, B-V, T e ff and log(g) on the one hand, and 
some abundance parameters, especially the one used to com¬ 
pute the tree, on the other. For instance, as a function of 
V, we see a decrease of [a/Fe] and possibly [Ti/Fe], and an 
increase of [Na/Fe], [Sc/Fe], [Lall/Fe]. For B-V, the most 
obvious trend is the increase of [a/Fe]. These correlations 
explain why OC1 is the only group to show a clear trend 
along its branch. 

The Na-0 diagram has been much scrutinized in the 
literature for identifying multiple stellar populations. It is 
shown in Fig[2] The Na-0 anticorrelation, which is present 
in to Centauri considered as a whole, is absent in individual 
groups, except in OC2a, and in the form of a slight gradient 
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Figure 4. Four projections of the parameter space in which the 
seven groups occupy distinct areas. 


in the OC1. [Na/Fe] is constant in OC3a and OC3b, and 
there is a correlation in OC3 c, w hich makes it identical to 
the third (metal-rich) group of GMttgr^gtjdJ (2Qlj|) and the 
fourth (metal-rich) group of ISimpson et al.l ( 2012ll . In the 
nomenclature of lCarretta et ah (120091') .' OC3 can roughly be 
associated with the extreme population, whereas the two 
other main groups would be a mixture of primordial and 
intermediate populations. 

Another noticeable anticorrelation among chemical ele¬ 
ments in w Centauri is that between [Ni/Fe] and [Fe/H], It is 
weak but significant (correlation coefficient -0.3, p-value=0), 
as seen on Fig. [3] Among individual groups it is present in 
OC3c and OC3b, and weakly in OC1. As mentioned ear¬ 
lier, OC2a and OC2b are only distinguished by [Ni/Fe]. 
The anticorrelation seems to be present wh en considering 
abundances integrated over the whole cluster. IGratton et all 
(l2004l'l show such a diagram, but dismiss as deviant the low 
[Ni/Fe] in their three metal-rich globular clusters. 

The seven groups occupy distinct locations on three pro¬ 
jections of the parameter space: [Na/Fe] vs [a/Fe], [Lall/Fe] 
vs [O/Fe] and La/Eu vs [Fe/H] (Fig. O. In the first pro¬ 
jection, OC1 has low [a/Fe], and there are three paral¬ 
lel diagonal sequences, OC2c, OC2a+OC2b+OC3a-|-OC3b 
and OC3c. In the second projection ([Lall/Fe] vs [O/Fe]), 
OC1, OC2a+OC2b and OC2c are of increasing [Lall/Fe] 
at high [O/Fe], OC3a and OC3b are of increasing [Lall/Fe] 
as [O/Fe] decreases, and OC3c has low [O/Fe] and high¬ 
est [Lall/Fe] . In the last projection, OC1 and OC2b have 
low La/Eu, OC2a shows a large spread, OC3a has interme¬ 
diate values of La/Eu, and the other groups (OC2c, OC3c 
and OC3b), which tend to be more metal-rich have a high 
La/Eu. 


Figure 5. Diagram of effective temperature vs gravity for the 
different groups. The colors of the groups are the same as in the 
previous figures. The full symbols correspond to [a/Fe] > 0.29. 
The lines are Dartmouth model isochrones for an age of 12.5Gyr. 
The metallicity and [a/Fe] values of the different isochrones are 
those indicated in Table [3] except in the top left box, where the 
limits at low metallicity are for [a/Fe] = 0.2 (solid line) and [a/Fe] 
= 0.0 (dotted line). The isochrones are all for a normal He abun¬ 
dance, except in the bottom left box (OC3b), where the limits at 
low metallicity are for Y = 0.245 (solid line) and Y = 0.4 (dotted 
line). The He abundances of the high metallicity isochrones for 
OC2c, OC3b and OC3c are obviously too low. 


3.2 Comparison with models : log(g)-T e fj diagram 

We now confront the chemical properties of the different 
groups with stellar evolution models. Our goal is essentially 
to determine whether Helium enrichment is necessary to ex¬ 
plain these properties, and, in that case, which groups are 
concerned. 

To t his end we use th e Dartmouth stellar evolution 
modelsQ dDotter et aklliooil ). which provide isochrones with 
different abundances of He and a-elements. The possible 
choices of He abundances are Y = 0.245+1.5Z, 0.33 and 
0.40. Unfortunately, the higher Y are all for [a/Fe] = 0. We 
plot these isochrones in the log(g)-T e g diagram, rather than 
in the usual color-magnitude diagram, because Teff is inde¬ 
pendent of distance. In view of the large spread in metallicity 
in most groups, we plot two isochrones corresponding to the 
two extreme values (i and f in Table [3J of [Fe/H] and [a/Fe], 
limits within which the data should in principle be included. 

We stress that the limits in [Fe/H] and [a/Fe] (given 
in Table [21 are not free parameters, and that the only free 
parameters are Y and age. We adopted an age of 12.5 Gyr 
for all the models; a different v alue wit hin the permissible 
range (11.4 to 13.1 Gyr, see Ijoo fe Leell2013l ) would not 
displace the isochrone significantly. 

The result of the model fitting is shown on Fig.0 OC1, 


1 http://stellar.dartmouth.edu/models/ 
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Table 3. The two extreme values (i and f) of the metallicity, 
[a/Fe] and Y for the seven groups as used in Fig. [5] Question 
marks indicate uncertain values, for lack of reliable isochrones. 


group 

[Fe/H] 

[a 

/Fe] 


Y 


i 

f 

i 

f 

i 

f 

OCl 

-2.00 

-1.50 

0.0 

0.4 

0.245 

0.247 

OC2a 

-2.00 

-1.50 

0.1 

0.4 

0.245 

0.247 

OC2b 

-2.00 

-1.50 

0.2 

0.4 

0.245 

0.247 

OC2c 

-1.75 

-1.00 

0.2 

0.6 

0.246 

0.33? 

OC3a 

-2.00 

-1.50 

0.2 

0.5 

0.245 

0.247 

OC3b 

-1.50 

-1.00 

0.2 

0.4 

0.40 

0.40? 

OC3c 

-1.75 

-0.50 

0.2 

0.6 

0.246 

0.40? 


OC2a, OC2b and OC3a are plotted together in the top 
left box, because their limits in metallicity are the same, 
and the differences in limiting [a/Fe] make little difference 
in the models, shown as solid ([a/Fe] = 0.2) and dotted 
([a/Fe] = 0.0) lines at the low metallicity limit. Most model 
isochrones are without He overabundances. We did plot an 
isochrone with Y = 0.4 (dotted line in the bottom left box 
of Fig. 02), which is an extrapolation from the isochrone 
with [a/Fe] = 0.0. For the high metallicity limits of OC2c, 
OC3b and OC3c, we decided that the limiting [a/Fe] are 
too high for obtaining a reliable extrapolation. We just note 
that raising Y lowers the surface gravity of the stars be¬ 
cause of the reduced electron scattering opacity, and thus 
moves the isochrone up in Fig. [5] at a given Teff. The con¬ 
sequence is that higher He abundances are required for the 
limit at high metallicity in groups OC2c, OC3b and OC3c. 
This confirm s previous identifications of He enhancement in 
this c luster (iNorr ij 20041: Piotto_etaL 20051 : iGratton et al.l 
l201ll : iKing et~ 20121 : Ijoo fe Led " 2013h . The latter au¬ 
thors predicted Y = 0.39-0.41 f or the three isochrones with 
[Fe/ H]> - 1.55. Grattoji et_al J (l201l|j used the data from 
I Johnson fe Pilachowski ( 20101 1 to estimate Y in their inter¬ 
mediate metallicity population, and found Y = 0.25 for their 
groups 5 and 1 (average [Fe/H]= -1.636), and Y = 0.347 for 
their groups 3 and 2a (average [Fe/H]= -1.416). 

This indicates that, at least in u> Centauri, Helium en¬ 
hancement occurs at metallicities higher than [Fe/H]= - 
1.5. This is in qualitative agreement with anal yses of other 
globu lar clusters: 47 Tuc ([Fe/H]= -0.72 dj_CnscienzoetJd] 
2010th NGC 2808 ([Fe/H]= -1.14 iD’Antona et al.l 120051: 


Piotto et _aL 200^: iMarino et al.ll2014th N GC 6388 and NGC 


6441 ([Fe/H]= -0.55 and -0.46; e.g. ICaloi fe D’Antonal 
l2007lh which are all very massive globular clus ters with likely 
He_enhancement, and of M3 ([Fe/H]= -1.50 fCatelan et al.l 
l2009h where no He enhancement has been found. This might 
of course be a coincidence, as so few globular clusters are 
available for comparison, and one cannot exclude that mass, 
for example, plays a larger role than metallicity in the He 
enhancement of stellar populations. There is also indirect 
evidence that the second generation stars, which seem to be 
present in all globular clusters, should be He-enhanced (see 
IGratton et al.l[20121 . and references therein). This limit of 
-1.50 should nevertheless be kept in mind until further no¬ 
tice. 








Figure 6. Density distribution of the radial distance (in arcmin) 
for the different groups. The densities are scaled to the density 
for OC1 which is shown in all boxes for comparison. 


3.3 Radial distribution 

The radial distribution of stars in uj Centauri has been 
studied by many authors, because it can give clues to the 
star formation history of the cluster. There is a consensus 
that the intermediate-metallicity RGB stars are more con¬ 
centrated than the metal-poor ones, but the most metal- 
rich ([Fe/H]> -0.6) RGB stars, also called RGB-a, are 


found to be as concentrated as the metal-intermediate ones 


by some authors 

Pancino et al. 2000; 1 Johnson et al. 20091; 

Bellini et al. 200E 

), and the least concentrated by others 

(Hilker & Richtler 

2000: Castellani et al. 20071. their Ta- 


ble 3). The situation is more nuanced in our groups, as 
shown on Fig [S] The least concentrated groups are OC2a 
and OC2c, followed by OC1 and OC2b. The OC3a, OC3b 
and OC3c groups are the most concentrated. The contradic¬ 
tion among previous results concerning the metal-rich pop¬ 
ulation is resolved by our clustering of this population into 
two groups (OC3b and OC3c) of different spatial distribu¬ 
tions. One explanation for the different spatial concentra¬ 
tion might be mass segregation, if the groups are dynami¬ 
cally relaxed, which may be the case if the stars have a suf¬ 
ficiently anisotropic velocity distribution (i.e. elongated or¬ 
bits). But the dynamical relaxation can be questioned and 
the mass range along the RGB covered by the sample is 
small. Following the suggestion of the referee, we propose 
rather that this group formed where it presently is located, 
and, for that same reason, is kinematically colder than the 
others._ _ 

[Pancino et alj (|200Clh have also found that the metal- 
poor population is elongated along the east-west (major 
axis) direction, whereas the two metal-r icher populations 
are el ongated in the opposite direction. iHilker fe Richtlen 
(l2000h note that the metal-richest stars are asymmetrically 
distributed around the center, with an excess toward the 
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south. We do find spatial asymmetries among the groups, 
but perhaps of marginal significance in view of the small 
numbers and statistical incompleteness of our sample. 


3.4 Kinematics 


The kinematics of the stars in the different groups are 
likely to provide clues about their formation. Furthermore, 
if the different groups have different kinematical proper¬ 
ties, our clustering results will be validated a posteriori. We 
thus collected kinematical data f or the stars of our sam ¬ 
ple. Such da t a are av ailable from ISuntzeff fe Kraftl Il996ll: 
iMavor et al.l (Il997fl: iReiins et al.1 (120061'): P ancino et all 
( 20071') : Sollima et al.l (1200911 : Johnson et al. ( 2008fh We did 
not use the kinematical data of van Loon et al.l d2007l J be¬ 
cause of their large quoted uncertainties (8 km/s), even 
though they would have added 11 stars to the kinemati¬ 
cal sample, three of which confirm the large mean rotation 
of group OC3b (see below). We obtained 1433 velocities and 
uncertainties for 650 stars. We computed the velocity zero- 
point of the diff erent samples, a nd then added 1 km/s to the 
velocitie s of IReii ns et al.l (I2006T). and subtracted 2 km/s t o 
those of lSuntzeff fe Kraftl ( 1996ll and lPancino et aD d2007h . 
Retaining only the stars that are in the 735-star sample and 
removing 14 velocities which showed a difference of 10 km/s 
or more with the mean of the others, we took the mean of the 
measurements for each star, weighted by their uncertainty. 
This produced a sample of 623 stars with velocities. 

We then adjusted a simple rotation curve by a weighted 
least-squares method to the radial velocities V r of each 
group, Vr — V 0 T 2 V m r m rsin(9 — Q 0 )l(r 2 m + r 2 ), where V a 
= 232.35 km/s is the mean radial velocity of the sample, 
r m is the radius at the maximum velocity of rotation V m , 
9 is the position angle measured from north through east, 
and 9 0 is the orientation of the axis of rotation with respect 
to the north-south direction. We also tried a Toomre rota¬ 
tion curve, which gave comparable results, but with slightly 
higher dispersions. We imposed an uncertainty of 0.5 km/s 
to stars with a single measurement and an uncertainty be¬ 
low that level because such low uncertainties, which strongly 
affect the fit, seem unrealistic given the spectral resolution. 

The kinematical parameters for the different groups are 
given in Table 2] N is the number of stars with kinematical 
data in the (sub)group, and er is the dispersion of velocities 
around the mean curve. The velocities and velocity disper¬ 
sions are in km/s, the radii in arcmin, the angles in degrees. 
The uncertainties on V m are all less than 0.1 km/s. For group 
OC3a, the search for the optimal parameters did not con¬ 
verge, because the rotation velocities in that group do not 
decrease near the origin. For that reason, we fit a rotation 
curve with constant velocity of rotation in concentric annuli, 
and, for that group, we list in Table [I] the mean velocity of 
rotation in different ranges of r (listed in col. 4, instead of 
r m ). For reference, we note that the mean peak rotation 
velocity of oj Centauri, determined with larger samples, is 8 
km/s, and is reached at about 8 arcmin dMerritt et al.lll997l : 
Ivan de Yen et ahll2006l 'l. 

The results show that the groups indeed have different 
kinematics. The groups OC1, OC2a and OC2c rotate like 
the bulk of the stars; OC3a, OC3b and OC2b rotate very 
fast in the inner region (where OC3a and OC3b are mostly 
located); OC3c shows very little rotation, if any. When sig- 


Table 4. The kinematical parameters of the seven groups, as¬ 
suming a variable velocity of rotation, except for OC3a 


Group 

N 

Vm 

E m 

do 

(7 

OCl 

194 

7.4 

11.4 

322 

8.06 

OC2a 

104 

7.1 

8.3 

354 

7.53 

OC2b 

38 

13.3 

4.0 

328 

8.93 

OC2c 

91 

8.2 

8.2 

347 

9.48 

OCl-t-OC2a+OC2c 

389 

7.4 

10.2 

337 

8.29 

OC3a 

8 

18.4 

r<3.5 

80 

13.3 

OC3a 

12 

16.6 

3.5^r<7 

12 

5.24 

OC3a 

10 

4.7 

7s?r 

306 

6.04 

OC3b 

46 

17.3 

2.0 

328 

8.49 

OC3c 

120 

3.6 

3.3 

6 

7.76 


nificant, rotation is roughly around the minor axis, which 
makes an accretion scenario unprobable, except perhaps for 
OC3a, which does not seem to be in dynamical equilibrium, 
as its rotation pattern is independent of the mass distribu¬ 
tion in the cluster. 

We are not the first to find that the su bpopula tions of c o 
Centauri have dif f erent kinematics. iNorris et alJ d 1997h and 
Ivan de Yen et ahl (l2006l f found that the metal-rich compo¬ 
nent, which is spatially more centrally concentrated, shows 
no evidence of systematic ro t ation. The se results were 
contra d icted by Pancino et al] (120071 1 and Ivan Loon et al.l 
( |2007h . iFerraro et ali (|2002l ) found that the metal-rich, so- 
called ”RGB-a” component shows a different proper-motion 
signature from the metal-poor and meta l-i nterme diat e com ¬ 
ponents. Thi s has b e en disputed bvlvan d e Yen e t ahl d2006h : 
iBellini et al.l (l2009l l ; lAnderson fc van der Marel ( 201Ch . In- 
cidentally, the latter found that the proper-motion disper¬ 
sion increases with decreasing stellar mass, suggesting that 
the core of ui Centauri is in the process of establishing energy 
equipartition. The explanation for these conflicting results 
is that one cannot segregate the different populations of t j 
Centauri on the basis of metallicity alone. 


4 DISCUSSION 

4.1 Primordial and later generations 

We now attempt to relate our groups to populations and/or 
successive generations of stars, and, to this end, look for 
evidence of external pollution in the chemical properties of 
the groups. 

The OC1 group is obviously a primordial population. It 
has the lowest mean value of most chemical elements, such 
as [Na/Fe], [Al/Fe], [Ni/Fe], [Lall/Fe], La/Eu, [a/Fe], and 
a high [O/Fe]. At variance with the other groups, there is 
no marked evolution of the abundances with rank within 
the group, except perhaps Ni and Eu. The main evolution 
with rank is a decrease in the V and K magnitudes and 
in T e ff and log(g), which were not used in the statistical 
analysis. The group has a spatial distribution that is close 
to the mean one in the cluster, and rotates as expected from 
the mass distribution of the cluster (Ivan de Ven et ahll2006l l. 
suggesting that it is dynamically relaxed. Finally, this group 
has 31% of the stars, which is about the expected percentage 
for the primordial population (see Sect. [I}. 

We find three hints of external pollution in the chemical 
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properties of the OC2 group: the presence of the Na-O an- 
ticorrelation and a large spread in [Lall/EuII] in the OC2a 
group, and a significantly larger [Ni/Fe] in the OC2b group. 
The Na-0 anticorrelation is generally attributed to pollution 
from H-burnin g at high temperature , in either fast rotating 
massive stars dDecressin et alJ 120071 ') or interm ediate-mass 
asym ptotic giant branch (AGB) stars dD’Antona fe Ventural 
l2007h . Ni is a heavier element than Fe, thus pro du ced in 
supernova nucleosynthesis. iTsuiimoto fc Shigevamal ( 2003 1 
attribute the increasing La/Eu in w Centauri to pollution 
by AGB stars after gas removal by supernovae of type II 
(hereafter SN II) at the end of a first episode of star for¬ 
mation. We note however that these two groups have the 
same O abundance as OC1, and that they are not markedly 
He-enriched, indicating that they formed from gas that was 
not much polluted. 

The OC2c group belongs to the OC2 group, but, like 
the OC3c group, differs markedly from the other groups of 
its main group, and these two outstanding groups should be 
discussed together. Both are significantly He-enriched (see 
Table [3j and have high [a/Fe]; this makes them likely to be¬ 
long to a secondary population. The difference is in [Al/Fe], 
low in OC2c and high in OC3c. Both groups are also spa¬ 
tially extended, while OC2a and OC3a are centrally concen¬ 
trated. 

The group OC3 is significantly Al-richer than the other 
groups. Al-rich stars must be second-generation stars be¬ 
cause oftheJughTemperatures necessary for the MgAl reac¬ 
tion (lLanger et al.l 1 19931 1 . which cannot be reached in low- 
mass stars. Furthermore, since He is the main outcome of hot 
H burning, these polluted stars must also be He-enriched. 
However, only OC3c and OC3b are significantly He-enriched 
(see Table [3]). The group OC3a, is very similar to OC2a 
and OC2b, except for [Al/Fe], which is the highest of all 
groups, and [O/Fe], which is low enough to make it part 
of an extreme secondary population. An Na-O correlation, 
rather than anticorrelation, is present in OC3c. It can be ex¬ 
plained by advocating super-AGB stars of masses between 
6.5 and 8 Mq, whose ejecta contribute to the formation of 
new stars about 120 t o 200 Myr after the first generation 
dD’Antona et al . ml. 

In summary, there is evidence for increasing external 
pollution in the stars of OC2a and OC3a, but the outstand¬ 
ing properties of OC2c and OC3c, as well as the properties 
of OC3a, intermediate between those of OC2a and OC3a, do 
not call for a simple division of the sample into three popu¬ 
lations, let alon e succ essive generations, of stars, as done by 
ICarretta et al.l (120091) for monometallic globular clusters. 


4.2 Origin of the enriched gas 

The next step is to use the spatial and kinematical proper¬ 
ties of the groups, as well as their chemical ones, to identify 
possible origins for the enriched gas that produced the sec¬ 
ondary populations. This gas was either expelled by winds 
from earlier stars or by supernovae. 

Gas ejected from stars with strong winds will disperse 
in interstellar space. Such stars can be either fast rotating 
massive stars or massive AGBs; the gas will thus be pro¬ 
duced a few tens or hundreds of Myr after the onset of the 
star formation event. Some of this gas may form stars when 
it is shocked during a passage of the cluster through the 


galactic plane (iOstriker et al.lfl972l '). before the gas density 
in the disc is high enough for ra m pressure stripping to occur 
(ITsuiimoto fc Shigevamall2003l'l . S i nce th e period of rotation 

1099), this plane 


of u> Centauri is 120 Mvr dDinescu et al. 


crossing is likely to happen every 60 Myr. iLeon et al.l (l200d 'l 
have found evidence for disc shocking in uj Centauri, in the 
form of a tidal tail aligned with the tidal field gradient. The 
stars formed in this way will show little mean rotation, if 
any. The remainder of the ejected gas will likely fall radi¬ 
ally into the strong potential well of the cluster, generating 
stars that will move in very elongated orbits. The gas that 
originates in stars with slower winds will presumably follow 
the stars in their orbits until close to the pericenter, where it 
collides with other gas clouds and settles into circular orbits, 
along which stars will form. 

This scenar io, which is inspired from the one proposed 
by lEggen et al.l dl962f) to interpret the properties of first- 
and second-generation dwarf stars in the Galaxy, can explain 
OC3. The group OC3c (chemically evolved and spatially ex¬ 
tended) shows little or no mean rotation, while OC3a and 
OC3b (chemically evolved and centrally concentrated) have 
the fastest m ean rot ation. This f aster rotation has also been 
predicted by iBekkil d20ld . l201ll l ; his hydrodynamical sim¬ 
ulations show that the gaseous ejecta from AGB stars pro¬ 
duce stars with different kinematics from the ones of the 
primordial population, namely a higher mean rotation and 
a lower mean velocity dispersion. In any case, OC3a and 
OC3b have to be fairly young, because this kinematical 
difference with the other groups cannot survive many Gi- 
gayears of dynamical friction. The core relaxa tion ti m e (at 
a radius of 2.5 aremin) of uj Centauri is 4 Gyr dHarri 3 119961 . 
2010 editior0), suggesting ages of that order. There are al¬ 
ternative explanations for faster rotation, either the gravo- 
gyr o instab ility, which affects massive stars more than others 
dKiiri et al]|2004h . or a central binary black-hole, which pro¬ 
duces ” suprathermal” stars at a few core radii dMapelli et alJ 
l2005h . but such mechanisms do not favor metal-rich stars 
over metal-poor ones. 

Supernova explosions are the most common source of 
enriched gas in the Galaxy. The progenitors of SN II are 
massive stars, whose lifetime is shorter than 100 Myr, and 
can be as short as a few Myr for the most massive ones. 
The progenitors of SN l a are binary stars; their peak of 
occur rence is 40-50 Myr (Matteucci fc Recchil[200l]j or 100 
Myr dKobavashi fc Nomotol 2009l l. However their frequency 
in globular clusters is probably much lower than that of SN 
II dVoss fe Nelemaiisll2012l : IWashabaugh fe Breemanll2013l l. 

Both types are more likely to occur in the central region of 
the cluster, where massive and binary stars are preferentially 
found. SN II produce O and a-elements, but little Fe, while 
SN la produce mainly 56 Ni (which decays to form 56 Fe), 
and other isotopes of Ni in lower proportions, b ut not a- 
elements, so that [a/Fe] actually decreases (|lwamoto et aid 
119991 . , their Table 3). SN la enrichment has been found in 
a few stars of the me tal-rich ”RGB- a” branch of uj Centaur i 
dPancino et al.lfioiil l. According to lMarcolini et ahl d2007t) . 
SN la occur in small pockets of gas before being mixed with 
the interstellar medium, at variance with SN II, which pol¬ 
lute the interstellar medium uniformly. Stars formed in these 


2 http://physwww.physics.mcmaster.ca/ harris/mwgc.dat 
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pockets are thus likely to retain the angular momentum of 
the progenitor gas clouds. This mechanism could explain the 
formation of OC2b, although its [a/Fe] is not lower than av¬ 
erage. 

The problem with enrichment by SN II is that these 
might be too energetic e vents for the gas to b e re¬ 
tained by the cluster (e .g. iTsuiimoto fc Shigevamal 120031 : 
iBaumgardt et all [2008 1. although radiative losses may 
prevent this from happening in a very massive cluster 
(|Marcolini et alJ 120071 ). SN events might not be an impor¬ 
tant source of enriched gas, thus perhaps not at the origin 
of the groups beyond OC1, but they can very well have trig¬ 
gered star formation by compressing and heating the gas 
enriched by other means. 


4.3 Formation scenario for the groups 

The formation scenario for the different groups of stars in 
oj Centauri must call on discrete events. This is required by 
the distinct isochrones found in the CMD of the cluster (see 
Sect, [lj, which lend credibility to our clustering of the stars 
into groups with distinct properties. Another requirement 
of the scenario is to take into account the large metallicity 
spread in the primordial and secondary populations, which is 
unique among galactic globular clusters. The current favored 
explanation is tha t the cluster is the remnant of a captured 
dwarf galaxy (e.g. iDinescu et al ] 19991 : iBekki fe Norrisll2006t 

iMarcolini et al] 120071 ). Our analysis has not enabled us to 
provide an explanation for the large metallicity spread in 
the primordial group. 

In the primordial population (OC1), the red-giant 
branch is the narrowest and the abundances are the least 
dispersed and the most constant with rank, compared to the 
other groups. These are properties of a population formed 
over a short period of time in a single event of star formation 
(a single generation, making a true isochrone) in a homoge¬ 
neous environment (very little spread in the evolutionary 
track). 

Judging from the chemical properties of the different 
groups and from their tracks on the CMD, OC2 and OC3 
must have formed successively from progressively enriched 
gas. We propose that OC2 formed from early SN explosions, 
before AGB pollution had time to enrich the gas, except per¬ 
haps for OC2a, in which there is some evidence for the Na-0 
anticorrelation. The group OC2b is enriched in O and Ni, 
which indicates pollution from SN la ejecta. This group is 
also more concentrated than OC2a and OC2c. In order to 
explain the high rotational velocity of this group, we have 
to assume that the SN la pockets formed near the center in 
or after the re-collapse phase at the end of the SN II events, 
thus acquiring kinetic energy. The group OC2c formed over 
the full extent of the cluster during early galactic plane cross¬ 
ings. 

The group OC3, which has the highest [Al/Fe] and 
[Na/Fe] and the lowest [O/Fe], formed last, in an interstel¬ 
lar medium O-depleted and Na-enriched, and also possibly 
He-enriched, by relatively massive AGB stars. The groups 
OC3a and O C3b formed by the mechanism proposed by 
iBekki! (i2010i . l201ll') in the central regions, with higher ro¬ 
tation velocities, while OC3c formed like OC2c. 

We realise that many questions are left unanswered by 


this scenario but it outlines possible reasons for the existence 
of seven distinct groups of stars. 


5 CONCLUSION 

We have performed a multivariate statistical analysis of the 
chemical properties of star samples in uj Centauri. At vari¬ 
ance with previous studies of the stellar populations of this 
cluster, we made no a priori assumptions about the num¬ 
ber of groups. In addition, our technique for obtaining the 
grouping is based on the relationships between the stars and 
not on their similarities, an approach that is suited to find 
populations and generations. Using the full range of known 
chemical properties, rather than just metal abundances, al¬ 
lows one to perform a clustering into groups that is much 
more efficient at revealing the distinctive properties of the 
stellar populations in a sample that at first sight shows a 
bewildering array of chemical, spatial and kinematical prop¬ 
erties. 

We found seven groups, which can be gathered into 
three main groups, and identified OC1 as the primordial 
population. The other groups show progressive chemical en¬ 
richment, while retaining distinct properties, for which we 
provide a tentative interpretation in terms of locus and 
mechanism of formation. Group OC2b formed near the cen¬ 
ter from SNIa ejecta and OC2a formed after AGB stars 
had chemically enriched the interstellar medium. OC3a and 
OC3b formed later, from the ejecta of stars with slow winds, 
OC3a being still dynamically unrelaxed. OC3c formed like 
OC2c, during Galactic plane crossings, but at a later stage 
of formation of the Galaxy. 

Our proposed clustering into groups lifts the contradic¬ 
tions between the results of earlier studies relative to the 
spatial distribution and kinematics of the different stellar 
populations, which cannot be distinguished by metallicity 
alone. It also proposes directions in which future investiga¬ 
tions could proceed to explain the uncovered properties. In 
particular, it remains to be shown whether the Na-O cor¬ 
relation found in OC3c can be produced by star formation 
during Galactic plane crossings, and whether the He-rich 
stars presently observed in globular clusters should on aver¬ 
age be more massive than if they had formed with a normal 
He abundance. 
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APPENDIX A: MAXIMUM PARSIMONY 
(CLADISTIC) ANALYSIS 

The cladistic analysis in astrophysics (astrocladis- 
tic) has been pre se nted in detai l in several p aper! ” 


(, FFaix^Burnet^e^al. 


2006al lbllcl: iFraix-Burnetl 12003 : 

iFraix-Burnet et~ l2009l l. We refer the reader to these 
references for a complete description of the method used 
here. 

Phylogenetic systematics relies on synapomorphies, 
which are derived character states, to infer common ancestry 
relationships. Such characters may be viewed as evolution¬ 
ary novelties appearing in a particular lineage. It is assumed 
that two closely relative objects share derived characters, 
which presumably originate in their common ancestor. 

An important advantage of phylogenetic systematics is 
to avoid the grouping of objects based on similarities due 
to evolutionary convergences or reversals. The success of a 
cladistic analysis depends on the behaviour of the parame¬ 
ters. In particular, it is sensitive to redundancies, incompati¬ 
bilities, and especially to homoplasies (reversals, parallel and 
convergent evolutions). The use of the parsimony principle 
minimises the number of homoplasies (see below). 

We have discarded Vmag and Teff since they are cer¬ 
tainly homoplasies (see Sect. not . The seven chemical pa¬ 
rameters used for the analysis can reasonably be consid¬ 
ered as synapomorphies since these abundances are gener¬ 
ated within stars before being ejected to form new stars. 
Hence they are specific to the ancestors that gave birth to 
a given population. Nevertheless, we cannot not guarantee 
them to be perfect synapomorphies because of measurement 
uncertainties and of our lack of a complete and thorough un¬ 
derstanding of the chemical evolution of the gas within stars 
and during ejection. 

The values for each parameter were discretized into 
30 equal-width bins representing supposedly evolutionary 


states. This choice of 30 bins is justified by a fair represen¬ 
tation of the diversity and of the continuous nature of the 
data, as well as a good stability of the analysis in the sense 
that the result does not depend on the number of bins. 

We imposed all the parameters to be ordered, i.e. 
changes between two adjacent states are more probable than 
between distant ones, independently from the sense of the 
change. We believe that this is a fair measure of the cost of 
evolution for continuous parameters. 

We adopted the popular parsimony criterion, which se¬ 
lects the most parsimonious tree among all possible arrange¬ 
ments because it represents the simplest evolutionary sce¬ 
nario compatible with the input data. This optimal tree 
corresponds to the minimum number of changes of states 
for all parameters that occur along the paths between all 
objects. This number is unique to each tree. 

The maximum parsimony searches were performed 
using the h euristic algori thm implemented in the 
PAUP*4.0bl0 dSwoffordll2003ll package. Heuristic methods 
do not explore the parameter space of all possible tree ar¬ 
rangements, which would take a prohibitive computer time 
for hundreds of objects, but try to find the minima. They 
cannot guarantee finding the absolute most parsimonious 
trees but generally require far less computer time while 
being quite effective. 

Fig [T] is a majority-rule consensus tree reflecting the 
most common features in all equally most parsimonious 
trees. 

The careful examination of the resulting tree must be 
performed to check a posteriori the behaviour of the param¬ 
eters. Firstly, the robustness of the tree is a good indicator 
of the quality of the characters (see below). Secondly, the 
evolution of the parameters along the tree can use for in¬ 
stance Fig. nm or Fig. IB2I Note that this should ideally be 
done with the binned parameters. Other analyses can then 
be repeated after removal of suspicious parameters. For in¬ 
stance, [Sc/Fe] and [Ni/Fe] seem to be too much variable 
on Fig. IB 1 1 An analysis without these two parameters still 
produce the three main groups. 

The robustness of the tree was assessed in the following 

way: 

(i) using several analyses modifying sligthly the set of pa¬ 
rameters. This is the philosophy of the bootstrap technique, 
which would be too much time consuming due to the large 
number of objects. In all cases, the result did not change 
much in the sense that the three main groups can be iden¬ 
tified quite easily. 

(ii) considering the quality of the consensus tree. A con¬ 
sensus is a summary of all the equally parsimonious trees 
found by the software. We found each node to be present at 
least in 89% of the 9517 equally parsimonious trees, most 
being in 99 or 100% cases. 


APPENDIX B: COMPLEMENTARY FIGURES 


3 see also http://astrocladistics.org 
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Figure B2. Boxplots 
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